Fast Relaxation in a Fragile Liquid under Pressure 
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The incoherent dynamic structure factor of ortho-terphenyl has been measured by neutron time- 
of-fiight and backscattering technique in the pressure range from 0.1 MPa to 240 MPa for tempera- 
tures between 301 K and 335 K. Tagged-particle correlations in the compressed liquid decay in two 
steps. The a-relaxation lineshape is independent of pressure, and the relaxation time proportional 
to viscosity. A kink in the amplitude /q(P) reveals the onset of (3 relaxation. The /3- relaxation 
regime can be described by the mode-coupling scaling function; amplitudes and time scales allow 
a consistent determination of the critical pressure P C (T). a and f3 relaxation depend in the same 
way on the thermodynamic state; close to the mode-coupling cross-over, this dependence can be 
parametrised by an effective coupling T oc nT~ 1//4 . 
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The glass transition can be induced by decreasing the 
temperature or by increasing the pressure. For practical 
reasons, most experimental investigations concentrate on 
temperature effects at ambient pressure Pq = 0.1 MPa, 
although there is a clear interest in studying supercooled 
liquid dynamics in the full P, T parameter space. Vari- 
able pressure measurements have been decisive in demon- 
strating that density is not the only driving force behind 
the glass transition jl]-[|. In the microscopic approach 
of mode-coupling theory 0] , both pressure and tempera- 
ture control the dynamics through variations of the static 
structure factor. Close to the dynamic cross-over, the ef- 
fect of these variations can be expressed through a single 
separation parameter which however is not given by den- 
sity alone, as was confirmed by depolarised light scatter- 
ing in the fragile liquid cumene |J. 

Here we report on incoherent neutron scattering in 
the fragile van der Waals liquid ortho-terphenyl (OTP, 
Tg(Po) = 243 K) at variable pressure. As function of 
temperature, the microscopic dynamics of OTP has been 
studied extensively by neutron scatter ing pi fj| , light 
scattering |||| and molecular dynamics [|l0|]ll[ |. As func- 
tion of pressure, viscosity fll2|l , photon correlation [jUg , 
dielectric loss [Q , and specific heat [||,[l5) data are avail- 
able for comparison. From the Debye- Waller facto r [p|JIo| 
and from direct observation of fast (3 relaxation | -|9| a 
cross-over was located at T C (P ) = 290 ± 5K @. 

OTP (Aldrich) was purified by repeated crystalli- 
sation from liquid methanol and vacuum destination. 
An A17049.T6 high pressure cell of the Institut Laue- 
Langevin (ILL) with a ratio of outer to inner diameter 
of 2 allowed us to attain up to 250 MPa at temperatures 
up to 335 K; helium was used to transmit the pressure. 

The experiments were performed using backscattering 
(IN16) and time-of-flight (IN6) spectrometers of the ILL. 
Incident wavelengths 6.27 A^ 1 and 5.12A -1 gave reso- 
lutions (fwhm) of about 1 ^eV and 80 ^eV, respectively. 



Vanadium was used to calibrate the detectors and to yield 
the resolution functions. The measured transmission of 
the empty pressure cell at Pq was about 84%, that of the 
sample about 88%. Self shielding and multiple scattering 
turned out be negligible at large scattering angles. This 
was demonstrated by comparison of spectra taken with 
and without the pressure cell at 290 K: the normalised 
intermediate scattering functions S(Q,t) fall onto each 
other. Therefore we base our conclusions exclusively on 
high Q data (> 1 A" 1 ). On IN16 we measured isoviscous 
and isothermal (320 K) spectra as well as the isothermal 
(301 K) elastic intensity. On IN6 we monitored two differ- 
ent isotherms (316 and 324 K). Preliminary results from 
isochoric measurements have been reported in [|l7| . 

Before turning to the fast dynamics, we address the 
long-time limit of a relaxation; results on a relaxation 
will be needed as input in the analysis of (3 relaxation. 

In a first survey, backscattering showed that the mean 
relaxation time follows the pressure dependence of viscos- 
ity. Then, the shape of a-relaxation spectra was deter- 
mined in a couple of long running scans. At three points 
(P, T) with equal viscosity, the quasielastic spectra all 
coincide (Fig. la) which confirms that the microscopic 
relaxation time is proportional to the macroscopic viscos- 
ity 7], and demonstrates that within experimental error 
the spectral distribution is independent of pressure. We 
deduce an average slope for isokinetic curves (dT/dP) T 
in the P,T plane of about 0.28±0.01 K/MPa. Similar 
values were found for the pressure dependence of the 
glass transition temperature (dT g /dP) jTa], and of the 
a-relaxation time in photon correlation p3| and dielec- 
tric spectroscopy Q. 

The invariance of the spectral shape is further illus- 
trated by construction of a master curve. The interme- 
diate scattering function S(Q, t) is obtained by Fourier 
deconvolution of the measured S{Q,uj). Since the a- 
relaxation time is proportional to viscosity, reduced times 
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can be defined via i = tr)(P )/r){P). The rescaled S(Q,t) 
coincide over two decades (Fig. lb). 

Time-temperature-pressure superposition is well es- 
tablished by phot on corr elation and dielectric loss spec- 

as well as in other fragile liq- 
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troscopy in OTP |13| 

uids jUI; for a given time shift, the spectral lineshape 
changes even less with pressure than with temperature 
jjjyiij. Our measurements extend the validity of this 
scaling principle to a picosecond range where contact can 
be made with fast molecular motion. 

Via the total intensity of the elastic or quasielastic line, 
backscattering may reveal the presence of faster relax- 
ation processes. In fact, the first manifestation of a dy- 
namic transition at T c was obtained from elastic scans 
j|. Here, we monitored the pressure dependence of the 
elastic scattering S(Q,uj ~ 0) at 301 K (Fig. 2). At 
high pressures, in the glassy phase, the elastic intensity 
is constant. This make the analysis even simpler than 
in temperature-dependent measurements where the har- 
monic evolution of the Debye- Waller factor had to be 
taken into account. On releasing the pressure, as on in- 
creasing the temperature, an anomalous decrease of the 
elastic intensity indicates the onset of j3 relaxation. This 
decrease is observed first and most pronounced for large 
wavenumbers. For the lowest pressures the elastic line 
begins to broaden, and its total intensity can no longer 
be determined from S(Q,uj ~ 0). Instead, one has to 
take the integral over the full a line. 

From the scaling analysis of isothermal a-relaxation 
spectra at 320 K (Fig. lb), we find that up to 105 MPa 
the integrated a-relaxation amplitude /q does not de- 
pend on pressure. This allows us to extrapolate /q at 
301 K from P towards higher pressures, depicted as hor- 
izontal lines in Fig. ||. The intersection with the elastic 
intensities (vertical line) which turns out to be indepen- 
dent of Q permits the identification of a critical pressure 
P c ~ 42 MPa. 

On IN6 we access the (3 process directly. As the 
quasielastic broadening is of the same order as the instru- 
ment's resolution width, data are analysed by Fourier de- 
convolution (Fig. 3). In our temperature-dependent ex- 
periments, spectra from different instruments were com- 
bined into an intermediate scattering function S(Q,t) 
which covered three decades in time ||[7j ; the fast j3 relax- 
ation could be described by the mode-coupling prediction 
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The dynamic range of the present study does not allow an 
independent verification of ([!]); nevertheless, taking ad- 
vantage of previous results, we can use (^) to extract the 
pressure dependence of Hq and t a from IN6 data. Be- 
sides the shape parameter A = 0.77 we input the Debye- 
Waller factor determined from a Kohlrausch fit to the 
long-time tail of low-pressure data. While the asymptotic 



law ([y) does not cover the short-time dynamics below 1 
psec, it consistently describes the bending of S(Q,t) into 
and out of the intermediate plateau /q (Fig. 3). 

Within mode-coupling theory, the a-relaxation ampli- 
tude /q, the /3-relaxation amplitude Hq, and the cross- 
over time t(j depend on the equilibrium structure fac- 
tor S(Q) which in turn is a regular function of para- 
meters like T and P. Close to the dynamic cross-over, 
this dependence can be expressed through a linear sep- 
aration parameter |po| , a oc (T — T c )/T c in isochoric or 
a oc (P c — P)/P c in isothermal experiments. On ap- 
proaching the cross-over from the liquid side, Jq varies 
only weakly, whereas Hq and t a are predicted to become 
singular: 
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with a critical exponent a = 0.295 determined by A. 

In Fig. 4a, b these power laws are tested for two 
isotherms by plotting Hq 2 and t~ 2a vs. P. Asymptot- 
ically, linear behaviour is observed, and extrapolation 
gives consistent estimates for P c at 316 and 324 K. Away 
from the asymptotic limit, the data scatter considerably, 
similarly as in temperature-dependent experiments . 

Together with T c at ambient pressure and the value 
P C (301K) obtained from elastic scans, a total of four 
points of the dynamic phase boundary T C {P) is ob- 
tained. These points, summarised in Tab. 1, fall onto 
a straight line with slope dT c /dP ~ 0.28K/MPa which 
is in remarkable accord with the slope determined from 
a-relaxation, viscosity and T g data: the mode-coupling 
cross-over line is parallel to lines of equal a response. 

Lines of constant density, on the other hand, have 
slopes of about 0.70K/MPa which confirms that 
temperature influences the dynamics not only via free 
volume. The same conclusion was reached by simulation 
of a Lennard- Jones system pl| and by light scattering in 
cumene ||; it was suggested instead that at least the a 
dynamics can be parametrised by viscosity [Bj . 

Our results imply that the dynamic cross-over can 
still be characterised by one single separation parameter 
tr(T, P) which however does not depend on density alone. 
Instead, it appears that density n and temperature T can 
be combined to an effective coupling V oc nT~ 1 / 4 : Tab. 1 
shows that T c is constant within experimental error. Fur- 
thermore, constraint fits to S(Q, t) confirm that the tem- 
perature and pressure dependence of amplitude and time 
scale can be described by the power laws (||) in conjunc- 
tion with just one separation parameter a = (T — T c )/T c . 

A coefficient T oc nT -1 / 4 is known to characterise equi- 
librium properties of a dense soft-sphere fluid with repul- 
sive core potentials of the Lennard- Jones r~ 12 type p2[ . 
In a supercooled liquid, we are actually in a high-density 
regime where the static structure factor S(Q) is sensi- 
tive only to the repulsive part of the potential 22 23] . In 
simulations of binary soft-sphere p4| and Lennard- Jones 
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p5[ systems, the glass transition was found to occur at 
constant T c . 

It was not expected that these results apply literally 
in a complex molecular liquid like OTP. In fact, the com- 
parison with a Lennard- Jones liquid can be made almost 
quantitative. To this end, we need the two parame- 
ters of the 6—12 potential: a van-der-Waals diameter 
a = 7.6 A, estimated from structural information [p6| , 
is in good accord with the recently observed prepeak in 
the static structure factor at Q ~ 0.85 . . .0.9 A" 1 |],|7|], 
and a potential depth e = 600 K has been determined 
by calibrating a three-site OTP model to experimental 
density and diffusivity fTofl , Using this input, we obtain 
r c = 1.50 (Tab. 1) which is closer to the Lennard-Jones 
result r LJ ~ 1.25 pl] , p5| than one could reasonably ex- 
pect. 

In conclusion, our experiments give a consistent pic- 
ture of microscopic dynamics in the T, P-planc. 

Time-temperature-pressure superposition, known from 
Hz-MHz spectroscopy, describes structural a relaxation 
down to the picosecond range where macroscopic prop- 
erties of the material evolve from its microscopic dynam- 
ics. Spectral lineshapes observed by incoherent neutron 
scattering are invariant, the mean relaxation time is pro- 
portional to the viscosity, and states of equal kinetic are 
connected by the same slope dT/dP, regardless of the 
time scale. 

Fast (3 relaxation is detected from the a-relaxation 
amplitude; its full spectral shape is observed by time- 
of-flight spectroscopy. Just as in variable temperature 
measurements, the data are fully compatible with mode- 
coupling predictions; they allow the unambigous and con- 
sistent determination of a dynamic cross-over line T C (P). 
This line is parallel to the lines of equal viscosity: a and 
P relaxation are driven by the same parameter. 

At least in a certain neighbourhood of T C (P) this pa- 
rameter is linear in the effective coupling T oc riT~ 1//4 
which in turn depends ultimately on the molecules' short- 
range interaction. 
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T(K) 


P c (MPa) 


n c (nm :i ) 


r c 


290 


0.1 


2.840 


1.495(4) 


301 


42(3) 


2.875 


1.499(6) 


316 


90(7) 


2.906 


1.497(5) 


324 


123(7) 


2.931 


1.501(9) 



TABLE I. Critical pressure P c and corresponding density n c J2i| as function of temperature. The product F c — na 3 (e/T) 1//4 
is found to be constant (its errors are estimated from AT C = ±3K at 290 K and AP C = ±5 MPa). The choice of a = 7.6 A and 
e = 600 K is explained in the text. 



FIG. 1. (a) Spectra from IN16 at Q=1.8 A -1 for three combinations of temperature and pressure which lead approximately 
to the same viscosity. The spectra are scaled to their value at u>=0. The line represents the measured resolution function. - 
(b) Master curve S(Q,t) from isothermal (320 K) data constructed by rescaling times with viscosity. The solid line is a fit with 
the Kohlrausch function /q exp(— t 13 ), with a stretching exponent (3 = 0.6. 

FIG. 2. Debye-Waller factor as a function of pressure at 301 K. Symbols show elastic scattering intensities: as long as there 
is no quasielastic broadening, they yield the Debye-Waller factor which is P independent in the glass (solid lines). Dashed 
horizontal lines indicate the integrated a intensity, measured at Po and shown to be P independent. The dotted vertical line 
indicates the critical pressure P c . 

FIG. 3. Intermediate scattering function S(Q,t) at Q=1.8A _1 from IN6 in the /3-relaxation regime for different pressures. 
Lines are fits with (^) with a fixed shape parameter A = 0.77 and fixed plateau values fg (indicated by thick circles). 

FIG. 4. /3-relaxation amplitude and cross-over time for two isotherms 316 and 324 K and three wave numbers Q=1.2 (□), 
1.5 (•), and 1.8 A -1 (o). The data are linearised as Hq and t CT 2a according to (H); lines are fits to Q averaged values. 



4 



reduced time (nsec) 
0.01 0.1 1 




Hco (|ieV) 



Tolle et al. Fig. 1 







~i n 1 r 



t 1 1 r 



~i r 



-0.25 



-0.5 



o 

o 

o ♦ 



o o o 



o 



o o 



o o 



♦ ♦ ♦ 





■ 






■ - 


■ 








-e— 


o 


u 




♦ 




♦ 


♦ 






♦ - 



O: 



O 



-0.75 



■ Q=1.2A" 1 
o Q=1.5A 1 
♦ Q=1.8A" 1 



o 



•1.25 



J L I L 



_l I I L_ 



_l L 







100 
P (MPa) 



200 



Tolle et al. Fig. 2 




Tolle et al. Fig. 3 




Tolle et al. Fig. 4a,b 



